:arXiv:astro-ph/0009399vl 25 Sep 2000| 



Table 4: List of all PNe with 8-13/im spectra, and the identified warm dust emission features. 



name 


I 


b 


dust 
type 


ref. 

a 


N/0 
type 


He 


c 


N 





ref. 

h 


D 

Z95 


D 

IRAS 


solar 












0.10 


3.6(-4) 


9.3(-5) 


7.4(-4) 


24 






SwStl 


1.50 


-6.70 





3 


III 


0.11 


3.9(-4) 


4.1(-5) 


5.9(-4) 


1 


3.79 


4.12 


NGC6369 


2.40 


5.80 


+ 


4 


I 


0.13 


5.8(-4) 


1.2(-4) 


2.9(-4) 


2 


0.75 


3.96 


Hb4 


3.10 


2.90 


+ 


7 


I 


0.13 


- 


2.8(-4) 


4.8(-4) 


2 


2.68 


9.52 


Hb6 


7.20 


1.80 


+ 


7 


I 


0.11 


- 


4.5(-4) 


5.1(-4) 


2 


2.45 


6.99 


NGC6309 


9.60 


14.80 


+ 


4 


II 


0.11 


1.3(-3) 


1.2(-4) 


9.5(-4) 


3 


2.35 


8.16 


NGC653r 


10.01 


00.7 


c 


- 


I 


0.19 


4.0(-5) 


8.9(-4) 


1.7(-4) 


23 


1.60 


3.41 


NGC6578 


10.80 


-1.80 


+ 


- 


II 


0.11 


- 


l.l(-4) 


5.6(-4) 


2 


2.47 


3.93 


M2-9 


10.80 


18.00 





3 


III 


0.08 


- 


5.0(-5) 


2.3(-4) 


4 


3.01 


2.82 


NGC6741 


33.80 


-2.60 


c 


7 


I 


0.11 


8.0 (-4) 


2.4(-4) 


5.4(-4) 


5 


2.33 


6.85 


NGC6572 


34.60 


11.80 


c 


3 


II 


0.09 


1.3(-4) 


9.6(-5) 


3.1(-4) 


6a 


1.02 


2.91 


NGC7009 


37.70 


-34.50 


+ 


4 


II 


0.10 


2.7(-4) 


1.8(-4) 


5.3(-4) 


7 


1.09 


4.35 


NGC6790 


37.80 


-6.30 


c 


1 


lib 


0.11 


7.2(-4) 


5.5(-5) 


7.6(-4) 


8 


3.18 


5.69 


CN3-1 


38.20 


12.00 


+ 


7 


III 


- 


5.0(-4) 


4.9(-5) 


1.3(-4) 


3 


4.05 


10.27 


NGC6210 


43.10 


37.70 


+ 


4 


lib 


0.10 


- 


4.9(-5) 


5.0(-4) 


3 


1.73 


6.84 


Vy2-2 


45.40 


-2.70 





3 


I 


0.09 


- 


4.4(-5) 


2.7(-4) 


4 


7.93 


3.93 


Hu2-1 


51.40 


9.60 


c 


7 


III 


0.09 


4.2(-4) 


4.0(-5) 


3.0(-4) 


3 


4.52 


10.50 


Ml-71 


55.50 


-0.50 


c 


- 


- 


- 


- 


- 


- 


- 


3.17 


5.90 


Hen2-447 


57.90 


-1.50 


c 


- 


- 


- 


- 


- 


- 


- 


6.30 


8.44 


IC4997 


58.30 


-10.90 





3 


Ill 


0.13 


4.1 (-4) 


2.0(-4) 


1.2(-3) 


6b 


4.43 


7.83 


NGC6886 


60.10 


-7.70 


C 


7 


II 


0.09 


4.1 (-4) 


1.3(-4) 


3.1(-4) 


9 


3.13 


10.67 


BD+303639 


64.70 


5.00 


c 


3 


lib 


- 


4.0(-4) 


l.l(-4) 


3.7(-4) 


10 


1.85 


2.12 


K3-53 


64.90 


-2.10 


c 


- 


- 


- 


- 


- 


- 


- 


6.95 


8.24 


K3-52 


67.90 


-0.20 


c 


- 


- 


- 


- 


- 


- 


- 


7.10 


7.51 


Hen2-459 


68.30 


-2.70 


C 


7 


- 


- 


- 


- 


- 


- 


6.02 


5.85 


M3-35 


71.60 


-2.30 


C 


- 


- 


- 


- 


- 


- 


- 


4.44 


6.85 


NGC6881 


74.50 


2.10 


C 


7 


I 


0.11 


- 


3.8(-4) 


5.5(-4) 


2 


3.96 


6.73 


NGC6884 


82.10 


7.00 


+ 


7 


lib 


0.10 


1.0(-3) 


3.0(-4) 


1.0(-3) 


11 


2.63 


9.23 


NGC6826 


83.50 


12.70 


+ 


4 


lib 


0.11 


1.7(-3) 


3.0(-5) 


2.4(-4) 


3 


1.42 


5.52 


NGC7027 


84.90 


-3.40 


c 


1 


II 


0.11 


6.9(-4) 


1.3(-4) 


3.1(-4) 


12 


0.64 


- 


Hul-2 


86.50 


-8.80 


+ 


- 


I 


0.16 


1.2(-4) 


3.2(-4) 


l.l(-4) 


3 


2.66 


17.49 


Ml-77 


89.30 


-2.20 


c 


- 


- 


- 


- 


~ 


- 


- 


5.29 


9.16 



comments 



silicatcs:0.95,203 UIR:0.05,2000 

weak continuum 

weak continuum 

weak continuum 

weak continuum 

graphite:0.20,157UIR:0. 56,0. 12,0. 12^ 

weak continuum'^ 

thick-silicates:0.38,536,1.09graphite:0.62,235 

graphite:0.29,128 UIR:0. 71,295 

SiC:0. 18,190 graphite:0.82,183 

weak continuum 

SiC:0.31,142 graphitc:0.69,188 

weak continuum 

weak continuum 

silicates:0.95,365 SiC:0.03,100 graphite:0.02,lll 

SiC:0. 17,200 graphite:0.83,175 

graphiteiO. 73,165 SiC:0. 27,186 

graphite:0.73,186 SiC:0. 27,179 

silicatcs:0.95,217 UIR:0.05,2000 

graphite:0.22,132 UIR:0. 78,355 

SiC:0.03,80 graphite:0.60,252 UIR:0.36,633 

graphite:0. 70,178 SiC:0.30,316 

graphite:0. 78,207 SiC:0. 22,262 

graphite:0.60,150 UIR:0.40,550 

graphitc:0. 24,191 UIR'':0.57,0.07,0.12 

graphite:0. 70,150 UIR:0.30,250 

weak continuum 

weak continuum 

SiC:0.08,151 graphitc:0.45,168 UIR:0.47,228 

weak continuum 

graphite:0.07,200 SiC:0.006,81 UIR'^:0.64,0.22,0.06 



Table 4: continued 



name 

IC5117 

21282S 

M2-49 

K3-62 

NGC6543 

K3-60 

Mc2-2 

IC5217 

B12-1 

M2-54 

K4-57 

Hbl2 

M2-56 

M4-18 

Ml-4 

IC2149 

K3-69 

Ml-5 

J900 

NGC2392 

Ml-6 

IC418 

IC2165 

Ml-11 

Ml-14 

Ml-12 

NGC3242 

IC2501 

Hcn2-47 



I 

89.80 

93.90 

95.10 

95.20 

96.40 

98.20 

100.00 

100.60 

104.10 

104.80 

107.40 

111.80 

118.40 

146.70 

147.40 

166.10 

170.70 

184.00 

194.20 

197.80 

211.20 

215.20 

221.30 

232.80 

234.90 

235.30 

261.00 

281.00 

285.60 



-5.10 
-0.10 
-2.00 

0.70 
29.90 

4.90 
-8.70 
-5.40 

1.00 
-6.70 
-0.60 
-2.80 

8.42 

7.60 
-2.30 
10.40 

4.60 
-2.10 

2.50 
17.30 
-3.50 
-24.20 
-12.30 
-4.70 
-1.40 
-3.90 
32.00 
-5.60 
-2.70 



dust 
type 

C 

C 

C 



C 

c 

+ 
C 

o 

+ 
o 
o 
c 



c 

c 

c 



c 

o 



rcf. 



C 



N 



O 



rcf. 

h 



N/0 He 

type 

II 0.11 7.9(-4) l.l(-4) 4.1(-4) 3 



lib 


0.11 


I 


0.13 


I 


0.14 



lib 



lib 
II 



lib 

II 

lib 

lib 

II 

III 

III 

III 

lib 

II 

lib 



0.10 
0.09 



0.12 

0.06 
0.07 
0.10 
0.02 
0.10 
0.02 
0.09 
0.11 
0.03 



1.3(-3) 



8.7(-5) 
1.9 (-4) 
6.9(-4) 



5. 6 (-4) 
4.3(-4) 
2.1(-4) 



1.3(-4) 



2.0(-3) 

2.1(-4) 

2.7(-4) 
4.1(-4) 



2.7(-4) 
6.1(-4) 



9.1 (-5) 
4.3(-5) 



4.3(-5) 
2.1 (-4) 
3.9(-5) 
8.3(-5) 
7.9(-5) 
2.7(-5) 
2.6(-5) 
3.2(-5) 
8.1(-5) 
1.6(-4) 
2.3(-4) 



3. 2 (-4) 
2. 9 (-4) 



3. 9 (-4) 
3.6(-4) 
3.4(-4) 
2.1 (-4) 
2.0(-4) 
5.5(-5) 
1.9(-4) 
1.3(-4) 
4.6(-4) 
5.1(-4) 
1.2(-3) 



0.10 l.l(-4) 6.0(-5) 2.2(-4) 13 8.11 



2 
14 



15 

16 

4 

6c 

17 

18 

19 

18 

3 

1 

4 



D D comments 

Z95 IRAS 

3.83 5.23 SiC:0.09,159 graphitc:0.61,213 UIR:0.30,369 

10.53 3.38 prominent UIR bands 

6.30 9.23 graphite:0. 12,330 SiC:0.01,70 UIR'^:0.34,0. 2,0.31 

4.07 8.98 graphite:0.79,224 SiC:0. 21,89 

1.12 3.36 weak continuum 

6.28 10.63 graphite:0.06,106 UIR'^iO. 59,0. 17,0. 18 

7.30 17.30 graphite:0.85,192 SiC:0.15,142 

- - weak continuum 

6.05 10.86 graphite:0.06,135 UIR'^:0.64,0.18,0.12 

13.80 11.45 silicates:1.0,173 

- - unidentified 
4.13 silicates:0. 71,304 SiC:0.23,2000 graphite:0.05,2000 
2.69 thick-silicates:0.53, 279, 1.32 graphite:0.47,2000 

6.85 10.53 SiC:0.03,232 graphite:0.43,137 UIR:0.54,113 

3.53 14.84 weak continuum 

2.07 8.85 unidentified (graphite at 157K) 

25.50 23.23 graphite:0.1 1,200 SiC:0.02,108 UIR'':0.48,0.21,0.18 

4.96 12.01 graphite:0.65,187 SiC:0.35,166 

3.40 11.04 graphite:0.34,143 UIR:0.66,133 

1.44 9.84 weak continuum 

4.35 10.50 SiC:0.21,110graphite:0.79,175 

1.02 2.62 SiC:0.26,120 graphite:0. 74,200 

2.52 11.94 SiC^:-,- graphite:-,90 

4.24 4.63 SiC:0.25,234 graphite:0.59,194 

4.37 13.57 weak continuum (possible SiC) 

6.48 10.11 SiC:0.24,140 graphite: 0.76,2 10 

0.94 5.51 weak continuum 

2.12 7.10 SiC:0.17,137graphite:0.83,202 

3.08 5.66 sihcates:0.66,242 graphite:0.34,143 



Table 4: continued 



IC2621 

NGC5315 

Hen2-131 

Hen2-117 

Hen2-113 

Hen2-142 

HDE330036 

Mz3 

CPD-568032 

Pcl-7 

NGC6302* 

Ml-26 

19w32 

Hb5 



291.60 
309.10 
315.10 
320.90 
321.50 
327.10 
330.80 
331.70 
332.92 
337.40 
349.50 
358.90 
359.20 
359.30 



-4.80 
-4.30 
-13.00 
2.20 
4.00 
-2.20 
4.10 
-1.00 
-9.91 
1.60 
1.00 
-0.70 
1.20 
-0.90 



dust 
type 

C 

C 

O 

-I- 

c 
c 
c 
o 
c 
c 
c 
o 
o 
c 



ref. 

a 

3 
3 
3 



N/0 
type 

I 
III 

lib 

I 
I 
lib 

I 
III 

I 



He 

0.09 
0.09 



0.18 



0.18 
0.08 

0.14 



C 

1.6(-3) 
1.4(-4) 

5.0(-3) 

1.9(-4) 

6.3(-3) 

1.0(-4) 
1.5(-4) 



N 



6.1(-4) 
1.8(-4) 



O 

8.3(-4) 
6.4(-4) 
4.7(-4) 



ref. 

h 

15 

15 
1 



6.6(-5) 4.8(-4) 22 



2.4(-4) 
3.5(-4) 
8.3(-5) 

8.3(-4) 
4.5(-5) 



2.3(-4) 
2.8(-4) 
4.8(-4) 

5.0(-4) 
3.0(-4) 



20 
16 
22 

21 
15 



l.l(-3) 6.6(-4) 



D 

Z95 
2.79 
1.96 
2.20 

5.50^' 
4.62 

1.17 

3.53^ 
0.57 
2.42 
4.30^ 
1.32 



D 

IRAS 
6.19 
4.15 
3.89 

1.95 
5.02 
5.05 
1.84 
1.85 
3.89 
1.59 
2.76 
4.38 
3.09 



comments 

SiC:0. 1,144 graphite:0.5,256 UIR:0.4,211 

SiC:0. 16,2000 graphiteiO. 55,271 UIR:0. 29,125 

silicates:0.42,409graphite:0.58,162 

weak continuum (possible UIR) 

UIR bands (like NGC 7027) 

graphite:0.11,140UIR'':0.64,0.16,0.09 

UIR bands 

thick-silicates:0.83,525, 1.42 graphite:0. 17,220 

UIR bands (like NGC 7027) 

graphite:0.18,152UIR'^:0.60,0.11,0.11 

graphite:0.66,154 UIR:0.34,490 

Silicates:l,201 

blackbody'':l,390,1.26 

graphite:0.41,162 UIR:0. 59,248 



" the references for the lO^m spectra are: 1) Aitken et al. 1979; 2) Aitken et al. 1980; 3) Aitken & Roche 1982; 4) Roche ct al. 1983b; 5) Roche & Aitken 

1983; 6) Roche et al. 1983a; 7) Roche & Aitken 1986; 8) Roche et al. 1991. 

'' Black body with optically thick foreground silicates with T(10/im)=1.26 

■^ free-free continuum subtracted, SiC and graphite are required but in a range of proportions and temperatures. 

"^ The UIR bands were fitted using individual resonance profiles. 

^ The spectrum of NGC 6578 presented here, with little continuum, is in contrast with that obtained with the IRAS LRS (Kwok et al. 1986), which shows 

prominent silicate emission. Observations at UKIRT show that a star ^30 " to the SE of NGC 6578 has a silicate emission band which has probably 

contaminated the IRAS LRS spectrum. 

^ This CGS3 spectrum was obtained with a 4" beam, which may account for the relative weakness of the 11.3/im feature compared to other nebulae (observed 

with 10" beams). 

B IRAS 21282+5050. The radio data for the S'(6cm)-based distance was taken from Likkel et al. (1994). 

'* The references for the gas phase abundances are: 1) Perinotto (1991); 2) AUer & Keyes (1987); 3) AUer & Czyzak (1983); 4) Koppen et al. (1991); 5) 

Hyung & AUer (1997); 6a) Hyung et al. (1994a); 6b) Hyung et al. (1994b); 6c) Hyung et al. (1994c); 7) Hyung & AUer (1995); 8) AUer et al. (1996); 9) 

Hyung et al. (1995); 10) AUer & Hyung (1995);11) Hyung et al. (1997); 12) Keyes et al. (1990); 13) Hyung & AUer (1996); 14) Feibelman et al. (1994); 15) 

Kingsburgh & Barlow (1994); 16) Zuckerman & AUer (1986); 17) Hyung (1994); 18) Cuisinier et al. (1996); 19) Costa et al. (1996); 20) Lutz (1984); 21) 

AUer et al. (1981); 22) De Marco et al. (1997).; 23) Feibelman et al. (1985); 24) Grevesse & Anders (1989) 

* The C abundances listed for NGC 6302 and NGC 6537 are probably underestimated: in NGC 6302, Cv is the most abundant ion, but was not included by 

AUer et al. (1981). We adopted C/O=0.9 for NGC 6302 (Mike Barlow, private communication), and did not use the C/0 information for NGC 6537. 

^ Radio data from the Strasbourg-ESO catalogue and references therein. 
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ABSTRACT 

We investigate the galactic disk distribution of a sample of planetary nebulae charac- 
terised in terms of their mid-infrared spectral features. The total number of galactic 
disk PNe with 8-13/Ltm spectra is brought up to 74 with the inclusion of 24 new ob- 
jects, whose spectra we present for the first time. 54 PNe have clearly identified warm 
dust emission features, and form a sample which we use to construct the distribu- 
tion of the C/0 chemical balance in galactic disk PNe. The dust emission features 
complement the information on the progenitor masses brought by the gas-phase N/0 
ratios: PNe with unidentified infrared emission bands have the highest N/0 ratios, 
while PNe with the silicate signature have either very high N enrichment or close to 
none, and SiC emission features coincide with a range of moderate N-enrichments. We 
find a trend for a decreasing proportion of 0-rich PNe towards the third and fourth 
galactic quadrants. Two independent distance scales confirmed that the proportion of 
0-rich PNe decreases from 30±9% inside the solar circle, to 14±7% outside. PNe with 
warm dust are also the youngest. PNe with no warm dust are uniformly distributed 
in C/0 and N/0 ratios, and do not appear to be confined to C/0~1. They also have 
higher 6 cm fluxes, as expected from more evolved PNe. We show that the IRAS fluxes 
are a good representation of the bolometric flux for warm-dust PNe. The requirement 
F(12//m) > 0.5 Jy should probe a good portion of the galactic disk, and the dominant 
selection effects are rooted in the PN catalogues. 
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1 INTRODUCTION 

Spectroscopy at 10 /im has brought significant informa- 
tion on the chemical composition of planetary nebu- 
lae (PNe). The dust signatures reflect the C/0 chem- 
ical balance at the tip of the asymptotic giant branch 
(AGB). In this article we use the 8-13/Ltm dust signa- 
tures as a systematic tool to investigate the distribu- 
tion of the C/0 abundance ratio in PNe. We classify 
the 8-13//m spectral signatures in a sample consisting 
of compact and IR-bright PNe in the Strasbourg-ESO 
catalogue (Acker et al. 1992). The sample excludes PNe 
traditionally associated with the galactic bulge, which 
will make the object of a forthcoming article. 

The family of emission bands usually referred to 
as the 'unidentified infrared bands' (UIR bands), with 
principal members at 3.3/im, 6.2^m, 7.7/im, 8.6/im and 



11.3/Ltm, were first observed in the mid-IR towards 
NGC7027 by Gillctt et al. (1973). Cohen ct al. (1986, 
1989) found good correlations between the strengths of 
all pairs of bands towards a sample of PNe, reflection 
nebulae and HII regions, showing that they correspond 
to a generic spectrum. The strength of the 7.7/im fea- 
ture relative to the total IRAS flux correlates strongly 
with the gas phase C/0 ratio in a sample of 6 PNe (Co- 
hen et al. 1986), and Duley & WiUiams (1981) identified 
the principal wavelengths of the UIR bands with tran- 
sitions in the chemical functional groups of aromatic 
molecules. Poly cyclic aromatic hydrocarbons are the 
commonly accepted carriers for the UIR bands (Leger 
& Puget 1984). Thus, although their exact carriers still 
remain to be determined, the UIR bands are indicative 
of a carbon rich environment. 

The 8-13/im spectra of PNe can also show signifi- 
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cant continuum emission attributed to 'warm' dust at 
~200K. As first shown by Aitken et al. (1979), a smooth 
emission feature with a peak at about 9.7/im is observed 
in the PNe SwSt 1, M 1-26 and Hb 12. The peak of emis- 
sion coincides with the Si-0 stretch in sihcates (Day 
1979, 1981), and this feature is similar to amorphous 
condensates of silicate materials (Day and Donn 1978). 
It is typical of the Trapezium region in Orion and of 
the circumstellar shells of some oxygen-rich stars (For- 
rest et al. 1975), and is also seen in absorption towards 
the BN infrared point source in Orion (Gillett et al. 
1975). Willner at al. (1979) detected a smooth emission 
feature with a rather flat profile, from 10.5/im to about 
12.7/im, towards the PNe IC418 and NGC6572. It is 
commonly observed towards C stars as excess emission 
over a black body spectrum, and is attributed to lattice 
vibrations in silicon carbide (Forrest et al. 1975). An- 
dersen et al. (1999) published the niid-IR transmission 
spectrum of meteoritic SiC grains; they obtained a good 
match to the C star feature when extracting grains with 
a small size distribution (i.e. < 5/xm). 

There are thus significant compositional differences 
in the dust content of PNe, which can be classified 
into dust emission types, and according to whether they 
contain 0-rich or C-rich grain materials (Aitken et al. 
1979). How do the lO/xm dust emission features compare 
with the gas phase abundances? What is the proportion 
of PNe that show each type of dust emission, and does 
their distribution show large scale variations across the 
galactic disk? It has been shown by Thronson et al. 
(1987) and Jura et al. (1989) that the proportion of 
C stars relative to M giants increases outside the solar 
circle, and it is interesting to investigate whether PNe 
follow a similar trend. In contrast to the C and M stars 
which span a range of locations in the giant branch, PN 
compositions reflect the surface abundances at the end 
of the AGB, with a well defined evolutionary status. 

This article is the first of a series devoted to the sta- 
tistical analysis of the grain composition in galactic disk 
PNe. In Section || we present a sample of 74 PNe with 
8-13/xm spectra, which includes 24 previously unpub- 
lished spectra obtained with CGS3 on UKIRT or the 
UCL spectrometer on UKIRT or the AAT. Section | 
also contains a brief description of the method used 
to classify the 10/^m continua (following Aitken et al. 
1979). We will then compare in Section the dust con- 
tent of PNe with their gas phase C/0 and N/0 abun- 
dance ratios, to show that the dust emission types rep- 
resent an alternative for determining the C/0 chemical 
balance and bring complementary information on the 
PN progenitors. The sky distribution of PN dust types 
is presented in Section 0. After adopting a statistical 
distance scale in Section ^, we will test in Section g 
the size of the PN sample in this work for stratification 
in Peimbert (1978) types, and discuss its homogeneity, 
giving some support for its statistical significance. The 
galactic disk distribution of the PN dust composition is 



Table 1: Log of observations. 





/ 


b 


instrument 


Date 




[dog 


:roes] 






NGC6537 


10.01 


0.70 


UKIRT+CGS3 


May 94 


NGC6578 


10.80 


-1.80 


UKIRT+UCLS 


Oct 87 


Ml-71 


55.50 


-0.50 


UKIRT+UCLS 


Jul 90 


Hen2-447 


57.90 


-1.50 


UKIRT+UCLS 


Jul 90 


K3-53 


64.90 


-2.10 


UKIRT+UCLS 


Jul 90 


K3-52 


67.90 


-0.20 


UKIRT+UCLS 


Jul 90 


M3-35 


71.60 


-2.30 


UKIRT+UCLS 


Jul 86 


Hul-2 


86.50 


-8.80 


UKIRT+UCLS 


Oct 87 


Ml-77 


89.30 


-2.20 


UKIRT+CGS3 


Sep 96 


M2-49 


95.10 


-2.00 


UKIRT+CGS3 


Sep 96 


K3-62 


95.20 


-1-0.70 


UKIRT+UCLS 


Jul 90 


K3-60 


98.20 


-1-4.90 


UKIRT+CGS3 


Sep 96 


Me2-2 


100.00 


-8.70 


UKIRT+UCLS 


Oct 87 


B12-1 


104.10 


-1-1.00 


UKIRT+CGS3 


Sep 96 


M2-54 


104.80 


-6.70 


UKIRT+CGS3 


Sep 96 


K4-57 


107.40 


-0.60 


UKIRT+CGS3 


Sep 96 


Ml-4 


147.40 


-2.30 


UKIRT+CGS3 


Sep 96 


IC2149 


166.10 


4-10.40 


UKIRT+CGS3 


Sep 96 


K3-69 


170.70 


+4.60 


UKIRT+CGS3 


Sep 96 


Ml-5 


184.00 


-2.10 


UKIRT+UCLS 


Oct 87 


Ml-14 


234.90 


-1.40 


UKIRT+CGS3 


Sep 96 


Hen2-117 


320.90 


-1-2.20 


AAT+UCLS 


Apr 87 


Hen2-142 


327.10 


-2.20 


AAT+UCLS 


Apr 86 


Pel-7 


337.40 


+1.60 


AAT+UCLS 


Apr 86 



presented in Section |7|. Section g summarises our con- 
clusions. 



2 8 13/im SPECTROSCOPY OF COMPACT 
AND INFRARED-BRIGHT PNE 

The criteria for selection of the PN sample were that 
they be compact, less than 10" in diameter so that 
most of the flux is contained in the spectrograph beam, 
and infrared bright, with IRAS 12/ini flux in excess of 
0.5 Jy. These criteria select the best candidates for the 
detection of the dust emission features. A sample of 
10 PNe was observed with CGS3 on UKIRT on 1996 
September 27 and 28, with the intention of increasing 
the number of objects measured beyond the solar circle. 
An observing log is presented in Table n^ where details of 
the previously unpublished spectra from the AAT and 
UKIRT are hsted. Both CGS3 and the UCLS were used 
in their low-resolution modes, producing oversampled 
spectra which are calibrated with respect to standard 
stars. Fluxes are accurate to 20%. The fluxes of emission 
lines present in some objects are listed in table 0. Al- 
though they provide information on the ionized gaseous 
component, we will not use the emission lines for the 
purpose of this work. The resulting spectra are shown 
in Figure |l|, together with fits to the 8-13/im continua 
based on the grain emissivities, which form the basis of 
the dust type classification. 

We classified the type of dust emission features 
according to the procedure described in Aitken et 
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Table 2: Emission line fluxes, in 10 ^'^W m 





[Arm] 


[Siv] 


[Ne ii] 




8.99 ^m 


10.52 /xm 


12.81 /im 


NGC6578 


2.10 


17.0 


— 


Ml-71 


3.57 


5.52 


2.27 


Hen2-447 


1.04 


— 


5.98 


K3-53 


— 


2.56 


— 


K3-52 


0.51 


2.45 


1.00 


M3-35 


— 


4.77 


1.99 


Hul-2 


— 


1.93 


— 


Ml-77 


— 


— 


0.98 


M2-49 


— 


2.29 


— 


K3-62 


1.16 


2.36 


0.58 


K3-60 


— 


2.27 


— 


Me2-2 


— 


0.63 


— 


M2-54 


— 


— 


0.90 


Ml-4 


— 


5.52 


— 


IC2149 


— 


— 


2.02 


Ml-14 


0.89 


— 


1.33 


Hen2-117 


4.03 


13.9 


3.92 


Hen2-142 


— 


— 


11.8 


Pel-7 


— 


— 


22.5 



al. (1979) and Aitken and Roche (1982). The emis- 
sivity functions e^ for the three types of grains are 
taken from the spectra of astrophysical sources, with 
F\ = eiB{\,T), where Fx is the observed flux density 
and B{X, T) is a Planck function at temperature T (e^ is 
then flxed by ei{10^in)—l). Additionally a smooth con- 
tinuum with emissivity /(A) oc A"-'^'*, which is taken to 
represent graphite grains (or amorphous carbon grains), 
was included in the fltting procedure. The relative con- 
tribution of emissivity functions and the temperatures 
of the grains are fit to the 8-13/ini continua following a 
X^ minimisation: 



^A = V aM>^)B{\,T,)/B{10^im,Ti), 

^ ^t 



(1) 

(2) 



where the sum extends to the number of dust compo- 
nents used in the fit, and ex stands for the error in 
each point of the observed spectrum _F'^'^^. The Planck 
functions are arbitrarily normalised at 10/im. Best fit 
values are thus obtained for the measure of contribu- 
tion of each emission type ai_4, and their black body 
temperatures Ti_4. The relative contribution of each 
dust type is a'^ = ai/J^i o-j, summing over all the dust 
types required in the fit. The coefficient a' is thus a rep- 
resentation of the fractional emission at A= 10/im. The 
fitting procedure assumes that all emitting materials 
are optically thin and that any absorption occurs in a 
cold foreground layer; in the case of absorption a factor 
g-'^W is included in Eq. g (t(A) is the wavelength de- 



Table 3: Classification of PN dust types. 

O (silicates) Silicate grain emission. 

C (UIR) UIR bands but not Sihcates. 

c (SiC) SiC grains, but neither Sihcates 

nor the UIR bands. 
+ (weak) weak or unidentifiable continuum. 



pendent opacity, with a profile given by the emissivity 
curve, keeping r(10/im) as a free parameter). The avail- 
able data provide insufficient constraints to warrant a 
more detailed radiative transfer treatment. 

As discussed in Aitken and Roche (1982), it is the 
silicate grain emission feature that really separates PNe 
into different groups. A portion of PNe with silicate 
emission also show the UIR bands, and sometimes re- 
quire SiC in small amounts. Graphite emission seems to 
be unrelated to the other types of grains. The PN dust 
signatures can be placed into four groups based on the 
dominant dust species at 10 /im as shown in Table l|. 

A list of all the identified dust emission features 
found in this sample of PNe can be found in Table 
^. The S/N ratio of many spectra are rather low, and 
higher quality data may confirm the need for mixture 
of grain types. However, in this article we are concerned 
with the dominant material. The column under 'com- 
ments' gives more information on the best fit parame- 
ters, in the form 



(graintype) : a',T,T, 



(3) 



where the optical depth field t is listed for only when 
absorption is required (which is the case for M 2-9, M 2- 
56 and 19w32). 

As would be expected from blackbody radiation be- 
tween 8-13//m, a typical dust temperature is >200K, 
and the 8-13/im dust emission can be referred to as 
'warm dust emission features' in comparison with colder 
dust <200K which makes the bulk of the FIR emission 
(e.g. Kwok et al. 1986). 

The objects are listed in the Strasbourg-ESO cata- 
logue (Acker et al. 1992) as 'true or probable' PNe. The 
exception is IRAS 21282-1-5050, whose optical spectrum 
(Cohen & Jones 1987) shows [O ill] emission lines, pho- 
tospheric absorption features corresponding to a heav- 
ily reddened [WC 11] nucleus, and a total luminosity of 
about 2x10"^ L0 for a guessed distance of 2 kpc, mak- 
ing it a probable low-excitation PN (the central star has 
been re-classified as an O star by Crowther et al. 1998). 

In the notation of Table S, out of 74 galactic disk 
PNe with 8-13/im spectra, there are 12 O nebulae (of 
which SwStl and IC4997 also show the UIR bands), 
16 c nebulae and 26 C nebulae (for 9 of which the fits 
are improved with the inclusion of SiC). The remainder 
shows either too little continuum emission (17 PNe), or 
no clear identification in terms of the classification used 
here (e.g. IC2149 and Ml-14 whose spectra are best fit- 
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Figure 1: 10/im spectra for the PNe listed in Table ||. In abscissae is the wavelength range in /im, and in ordinates 
the flux density in 10^^^ Wcm^^ /im"'^. A fit to the continuum emission is shown for the spectra with identified 
dust features, with parameters listed in Table ^. The data points around bright emission lines ([Am], [Siv] and 
[Neil] at 8.99, 10.52 and 12.81 /im respectively) are excluded from the fits. 
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ted by graphite emisssion only, and K4-57, whose spec- 
trum is flat in units of Janskys, and is atypical of PNe). 
Finally the fits in Hbl2 and Vy2-2, both classified as 
'O' PNe, require some amount of SiC, which is prob- 
ably due to variations in the silicate emission profile 
rather than a superposition of grain types. 



3 COMPARISON OF THE DUST 

COMPOSITIONS AND THE GAS PHASE 
ABUNDANCES 

How does the warm dust emission feature classification, 
based on the grain C/0 chemical balance, compare with 
the gas phase C/0 abundance ratios? How does it com- 
pare to nitrogen enrichment? In order to address these 
questions, we searched the literature for published gas 
phase abundances in PNe with 8-13/^m spectra. Table 
^ lists the PNe for which a detailed spectroscopic anal- 
ysis is available, but it should be borne in mind that 
the uncertainties in the abundance analysis are often 
substantial. The PNe are classified according to the Pe- 
imbert (1978) types, with the sub-types in N/O ratio 
introduced by Faiindez-Abans & Maciel (1987). Many 
assignments to N/O types are taken from the compila- 
tion in Maciel k Dutra (1992). 

Figure shows the distribution of gas phase C/O 
ratios for each 8-13/im continuum class. The correspon- 
dence is not direct, some nebulae with C based grains 
have gas phase C/0<1. There is however a good cor- 
relation in C based grains with C/O ratio, previously 
obtained by Barlow (1983) and Roche (1989), for ex- 
ample, but with a lower number of PNe. Silicate grains 
indicate an O rich environment, SiC a C rich environ- 
ment, and the UIR bands correlate with a strong over- 
abundance of C relative to O. Weak continuum PNe are 
widely spread in C/O ratios, suggesting they may corre- 
spond to later evolutionary stages, rather than C/0^1. 
It is established that the dust grain composition re- 
flects on average the gas phase composition, and can 
be used as probes of C enrichment in PNe . The warm 
dust features are thus an alternative to the UV lines 
[Cii] A2326, [Cm] A1908, [Civ] A1550, for coarse classi- 
flcations of C/0 ratios. 

The existence of 5 PNe with C/0<1 and grain emis- 
sion characteristic of C rich environments merits further 
attention. The same apparently paradoxical situation is 
hinted at by the superposition of silicates and C based 
grains in the spectra of SwStl, IC4997, Hbl2 (confirmed 
by measurements at 3 /im showing the 3.3 /im UIR 
band, Roche et al. 1996). Standard equilibrium chem- 
istry suggests the least abundant of C or O would be 
locked in CO (Gilman 1969), and the production of 
C-rich molecules in 0-rich circumstellar environments 
may be a rather rare non-equilibrium phenomenon (ob- 
served towards certain M supergiants in h and x P^r, 
Sylvester et al. 1998). This is also suggested by the ap- 
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Figure 2: Histogram showing the number of nebulae of 
each dust types as a function of gas phase C/0 ratios. 
Typical error bars are ~ -l-0.il/ — 0.15 horizontally for a 
20 per cent uncertainty in the abundance calculations. 

pearance of both C- and 0-type grains in novae (e.g. 
Smith et a. 1995). Thus mixtures of grain types imply 
either a stratification in the ejecta composition of the 
progenitor star, or the mixing of the progenitor ejecta 
with pre-existing material. This point is relevant when 
linking the grain types with the C/0 ratios of the pro- 
genitor stars, and will be discussed further in a forth- 
coming article (Casassus & Roche 2000, paper II). 

The link between the dust types and N enrichment 
is shown in Figure ||. The UIR bands are found mostly 
in nebulae of Peimbert type I, whereas SiC emission and 
'weak' continuum PNe are uniformly represented. Sili- 
cates are found either with strong N enrichment (type 
I), or none at all (types lib and HI). These trends should 
be confirmed with a larger sample of objects with known 
N/O ratios, especially for 0-type signatures. 

The stratification of PNe in N/O ratios with height 
above the galactic plane suggests that the Peimbert 
classification is indicative of progenitor mass. Thus, on a 
relative mass scale, the UIR bands correspond to higher 
progenitor masses, SiC to intermediate masses, silicates 
are found mainly for low mass progenitors, but also 
for the most massive ones. The uniform distribution of 
'weak' continuum PNe again suggests they may corre- 
spond to later evolutionary stages. 

But the dichotomy between O- and C-type PNe for 
high gas phase N/O abundance ratios indicates that 
there is no simple correspondence between progenitor 
mass and dust signature. The dust emission features 
provide complementary information to the Peimbert 
types. 



The Galactic Disk Distribution of Planetary Nebulae With Warm Dust Emission Features: I 7 



UIR 



Sic silicates weak 



1 — 1 




— 




— 


1 — 1 




■ 











No lib III I No lib II 



No lib III I No lib 



Figure 3: Relative proportion of PN dust types among 
Peimbert (1978) types, with sub-types from Faundez- 
Abans & Maciel (1987). 
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Figure 4: Sky distribution of PN dust types. 

4 THE SKY DISTRIBUTION OF PNE 
DUST TYPES 

The sky distribution of PN dust types, built using Table 
^, is shown in Figure |4[ Table [^ contains the moments of 
the distribution. We assumed the main source of errors 
on the fraction of each dust type is due to 'counting 
noise', and the uncertainties in the mean, < b >, and 
spread, 5rms, are estimated under the assumption that 
the parent distribution of PN latitudes is normal. 

The mean b is reasonably close to the galactic plane 
for all identified dust types. The much larger spread in 
'weak' PNe suggests they are closer on average; a larger 
spread in height above the plane is discarded on the 
basis that 'weak' PNe are uniformly distributed in Pe- 
imbert types, although the adoption of a distance scale 
is required to further discuss this issue. It seems that 
UIR PNe are the most concentrated towards the galac- 
tic plane, although SiC and silicate PNe have roughly 
the same spread in galactic latitude. For the subset of 
PNe with C based grains (the SiC and UIR nebulae), 
^rms=6.0±0.7, giving the ratio of the spreads for O rich 
PNe to C rich PNe, 6rms(0)/fe„ns(C) = 1.3 ± 0.3. AU 
uncertainties quoted in this article are ±lcr. 



Table 5: The properties of the sky distribution of plan- 
etary nebulae dust types. The errors quoted correspond 
to one standard deviation. 





<6> 


f^iins 


N 


l<90 


90 <; 




[degrees] 


[degrees] 




l>270 


l< 180 





-1.6±2.3 


8. Oil. 6 


12 


9 


3 


c 


-3.3±2.0 


7.9±1.4 


16 


8 


8 


C 


-0.4±0.8 


4.1±0.6 


26 


19 


7 


+ 


6.6±3.5 


15.5±2.5 


20 


12 


8 



There is a hint of a decrease in the proportion of 
silicate grains PNe from the first and fourth quadrants 
(—90 < I < 90) to the second and third quadrants 
(90 < Z < 270), from 0.25±0.07 to 0.17±0.09. Since 
it corresponds to only one a it cannot be considered a 
solid property of the distribution. There is also an in- 
dication of an increase in the relative proportion of SiC 
PNe, which doubles from 0.22±0.07 to 0.44±0.11. 



5 ADOPTED DISTANCE SCALES 

5.1 Distance scales based on 6 cm continuum 
emission 

One of the persistent problems related to the study of 
PNe is the difhculty of obtaining accurate distances. A 
review of the methods based on individual properties of 
PNe can be found in Peimbert (1992). These so-called 
'direct' distance estimators are available for a restricted 
number of objects, and suffer from intrinsic uncertain- 
ties. However, the distances to PNe can be determined 
on a statistical basis, which match the average proper- 
ties of PNe. A distance scale for PNe can be built on 
the assumption that a general relationship holds for a 
set of PNe. We adopted the distance scale derived by 
Zhang (1995), who used an arithmetic average of two 
complementary methods, one based on the mass-radius 
relationship (e.g. the review by Kwok 1994), and the 
other on the relationship between T^, the free- free 6 cm 
brightness temperature, and nebular radius (as intro- 
duced by Van de Steene & Zijlstra 1994). Zhang (1995) 
calibrated the mass-radius and Th-radius relationships 
with a large sample of PNe with individually and 'di- 
rectly' determined distances. This 'direct' method is ex- 
plained in detail in Zhang and Kwok (1993) and Zhang 
(1993), and depends on the distance-independent pa- 
rameters Tfc and the central star temperature T^,. Dis- 
tances thus obtained are strongly model dependent, and 
can be in disagreement with the more accurate com- 
parison of angular expansion rate and radial velocity. 
NGC 6572, NGC 6302, NGC 3242, NGC 2392, and 
NGC 7662 are given 'direct' distances of, respectively, 
2.9 kpc, 0.1 kpc, 1.1 kpc, 0.5 kpc, 1.6 kpc, while their ex- 
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pansion distance is 1.5±0.5, 1.6±0.6 kpc, 0.4±0.1 kpc, 
> 1.4 kpc, 0.8 ± 0.7 kpc (Gomez et al. 1993, Hajian et 
al. 1995, Hajian & Terzian 1996). But the details of the 
distances to each nebula is of secondary importance as 
long as the global properties of PNe are reproduced. 
In that sense, the Zhang (1995) distance scale gives 
a Gaussian distribution about the galactic centre for 
bulge PNe, with a narrower scatter than the scale by 
Van de Steene & Zijlstra (1994). 
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5.2 Distances from IRAS fluxes to optically 
thick PNe 

As the central stars of PNe evolve rapidly in time, and 
the nebulae are active radiatively and dynamically (e.g. 
Kwok 1994), distance scales based on invariant proper- 
ties of PNe cannot be applied to the whole PN popula- 
tion. But in the case of the sample discussed here, the 
4 IRAS band fluxes, coupled with constant luminosity, 
may provide an alternative distance scale, as we now 
argue. 

The compact and IR-bright PNe are likely to be 
young, surrounded by substantial molecular material 
and therefore optically thick to the ionizing radiation 
from the central star. In this case the total luminosity 
of the nuclei can be inferred from the flux of any HI re- 
combination line, by equating the number of H+ recom- 
binations to the number of photoionizations. Mendez 
et al. (1992) tested this hypothesis by comparing with 
spectroscopic studies of PN nuclei, linking the sur- 
face gravity and effective temperature to the luminosity 
through atmosphere models. Their conclusion is that 
most PNe are optically thin. However, their sample is 
biased against obscured central stars: out of 23 PNe, 6 
are infrared-bright and are among the sample discussed 
here, 4 of which have no warm dust. Thus Mendez et al. 
(1992) included only two nebulae with 8-13/im spectra 
showing warm dust emission, for which the ratio of lu- 
minosities derived from optical thickness to the model- 
atmosphere luminosities are 0.70 and 0.92 (for Ml-26 
and IC418). It is thus likely that the compact and IR- 
bright PNe with warm dust emission are optically thick. 

In PNe which are optically thick in the Lyman 
continuum, the ionized central regions are surrounded 
by substantial amounts of neutral gas, an environment 
favourable to dust-grain survival. Most of the UV radia- 
tion escaping from the ionized region would be absorbed 
by dust grains, which heat up as a result to ^100-200 K, 
and re-radiate in the mid- and far-IR spectral range. 
The IRAS band fluxes should give a good representa- 
tion of the bolometric fluxes using 

4 

where the sum extends to the 4 IRAS bands. 

PNe initially evolve at a constant luminosity once 



Figure 5: a) Synthetic PN luminosity function from 
the initial-final mass relationship of Groenewegen & 
de Jong (1993), and with progenitor ZAMS masses 
between 1.2 and 7 M©. b) The relationship between 
DiRAS and Zhang (1995) for the PNe with warm dust 
emission. 



they leave the AGB, as was first shown by Paczyhski 
(1970, see also Blocker 1995). The luminosity function 
for the youngest PNe should be close to that of tip-of- 
the-AGB objects. The distribution of core masses for 
stars at the tip of the AGB can be calculated with a 
synthetic AGB model and a crude galactic disk model 
(we used the analytic prescriptions in Groenewegen & 
de Jong 1993, and a galactic disk model described in 
paper II). The core-mass luminosity relationship from 
Wagenhuber & Groenewegen (1999), in the case of post- 
AGB objects (i.e. in the asymptotic regime and van- 
ishing envelope mass), gives the luminosity function 
of young PNeQ shown in Figure gja, with a mean of 
8500 Lq. a very similar luminosity function, with an av- 
erage of 9300 Lq , is obtained using the prescriptions in 
Wagenhuber & Groenewegen (1999) for the initial-final 
mass and core-mass-luminosity relations, and taking so- 
lar metallicity and an IMF index of 1 (instead of 1.72 in 
paper II, in a notation where the Salpeter (1955) IMF 
would be 1.35), with a constant star formation rate. 

It appears the PN luminosities are not expected to 
vary over more than one order of magnitude. Assign- 
ing the same luminosity of 8500 Lq for all PNe gives 
a maximum error on the distance of only a factor < 2. 
Distances to compact and IR bright PNe can thus be 
estimated under the assumption of constant luminosity. 
Eq. Q for the bolometric flux is likely to be a lower limit 
only, but in this article we are interested in the relative 
properties of the distribution of the different PNe dust 
types, and their absolute distances arc not required. We 
will refer to distances derived in this way by -D/_ras, and 
those from Zhang as Dzgs- 

We stress Djbas distances are only meant to in- 
vestigate an independent distance scale and its conse- 



* PNe progenitors for the sample discussed here were assumed 
to have masses in the range 1.2<M/Mq <7, see paper II. 
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quences on the derived galactocentric trends; these dis- 
tances should not be taken as accurate. 



5.3 Adopted distances to compact and 
IR-bright PNe 

The distances derived from the two methods described 
above are Usted in Table 0. Djuas distances appear to 
be reasonable for PNe with detected warm dust emis- 
sion. Figure ^ shows a good correlation between the 
two distance estimates in the case of PNe with warm 
dust emission: The ratio DirasI Dz95 is 1.87 on av- 
erage, with a l-CT spread of 0.84. This suggests that 
the luminosity used to derive Djuas distances may be 
overestimated by a factor '^ 3 — 4, if Dz95 distances 
are rehable. Also, the cases of NGC 6302, NGC 6572 
and BD-l-303639 (three PNe with warm dust emission 
features) allow comparing their expansion distances of 
1.6±0.6 kpc, 1.5±0.5 kpc and 2.68±0.81 kpc (Gomez et 
al. 1993, Hajian et al. 1995), with their Djjias distances 
of 1.6 kpc, 2.9 kpc and 2.12 kpc. Although the com- 
parison supports DjuAS distances, a handful of objects 
does not permit a generalization. In any case, DjfiAS 
may be used as an upper limit, except for PNe with 
upper limits in the IRAS 100/im band. It is worth not- 
ing, however, that the PNe Vy2-2, IRAS21282-f5050, 
Hbl2 and Hen2-113 are given distances on the Zhang 
(1995) scale that are in excess of Diras by a factor 
larger than 1.7 (which takes into account the maximum 
range expected in the PN luminosity function). 

In the remainder of this article we assign Djjias 
distances to PNe without radio data (i.e. M2-56 and 
HDE330036). The case of K3-69 seems to be anoma- 
lous: both distance estimates give ~25kpc, putting K3- 
69 at 1.7 kpc above the galactic plane, and we preferred 
to use the distance of 7.9 kpc from Cahn et al. (1992). 
Another anomalous case is M2-54: again, the Zhang 
(1995) distance scale and Diras both give a distance of 
~13kpc, placing it near the northern galactic warp. In 
order to avoid the uncertainties associated with exag- 
geratedly large distances, we adopted a maximum galac- 
tocentric radius of 14 kpc to compute the moments of 
the vertical distribution, thus excluding K 3-69 and M 2- 
54. 



6 TESTS FOR THE COMPLETENESS AND 
HOMOGENEITY OF THE COMPACT 
AND IR-BRIGHT PN SAMPLE 

The sequence in Peimbert (1978) types is indicative 
of progenitor mass, the highest being associated with 
type I. There is a stratification in height above the 
galactic plane as a function of N/0 type (e.g. Maciel 
and Dutra 1992). Such a stratification is indeed present 
in this sample. Inside the solar circle, the root mean 
square height over the plane, Zrms, is 0.14 kpc for type 



I, 0.37 for type Ila, 0.45 for type lib, and 0.68 for type 
III. This stratification is an indication that a statistical 
study based on the compact and IR-bright PN sam- 
ple would be sensitive to PN properties with the same 
dependence on progenitor mass as the Peimbert types. 
Out of 49 compact and IR-bright PNe with known N/O 
ratio, 29±6% are type I, 21±6% are type Ila, which is 
typical of PN catalogues (e.g. Maciel & Dutra 1992). 
The disk PN population seems to be homogeneously 
sampled, although the constraints will remain loose un- 
til a larger sample is available. 

A discussion of the selection effects is possible in 
terms of a comparison between the fraction of O-rich 
PNe and the predictions of synthetic AGB models. In 
paper II we compare expected tip-of-the-AGB statistics 
for the C/0 chemical balance with those from the dust 
signatures. We find that for a minimum PN progenitor 
mass of 1 Mq , about 50% of all young PNe should be 
O-rich, whereas we report 22%. To match the observed 
ratio, the minimum PN progenitor mass for the sam- 
ple in this work must be at least Af"""=1.2, at a 2a 
confidence level - assuming the warm dust composition 
correponds to the last ~ 2000 yr of AGB evolution. Av- 
eraging the AGB ejecta over the last 25000 yr increases 
the fraction of O-rich PNe by ~ 10%, which may explain 
the higher frequency of O-rich nebulae with plasma di- 
agnostics (40 ± 8% in the sample used here) . 

It is possible, however, that differences in the opac- 
ity functions among C- or 0-based grains could lead to 
different lifetimes of the warm-dust emission phase. In 
this case the tip-of-the-AGB and observed C/0 statis- 
tics would be different. To summarise, the maximum 
separation from the central star, ro, required to keep a 
dust grain at a temperature T > Tq is rp oc ^Jk^Jkr, 
where fc^ is the opacity averaged over the central star 
spectrum, and kx is averaged over a black body at Tq. 
For instance, it may be thought that if C-rich grains 
have higher k^^/kx than O-rich grains, then the C-rich 
warm-dust phase would be longer. But either A) kx is 
fixed, and then the acceleration of C-rich grains by UV- 
radiation pressure would be higher, thereby shortening 
the C-rich warm-dust phase, or B) fc^, is fixed and kx is 
lower for C-rich grains, in which case the lOfim fluxes 
of C-rich PNe would be lower, limiting the number of 
C-rich PNe in spite of their hypothetical extended life- 
times. In addition the role of the stellar wind and the 
interaction with the nebula considerably complicate the 
picture. 

But a test can be found for the preferential selection 
of O- or C-rich grains. The good agreement between 
DiRAS a-nd -Dz95 suggests Eq. Hcan be used to estimate 
the fraction of luminosity radiated in the 12/xm IRAS 
band by PNe with warm dust emission, with 



L 



12fim 



iyli'{12fim) 



(5) 



Omitting the PNe with upper limits in the 100/im IRAS 
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band, the fraction of luminosity emitted at 12/im is 
-0.25±0.15 for PNe with sihcate emission, 0.22±0.09 
for SiC PNe, 0.27±0.14 for UIR PNe, and ~0.25±0.14 
for all warm dust types, while 'weak' PNe have a 
Li2fim/L-k ratio of 0.11±0.03 (this sample is biased to- 
wards high values of Li2fj,m/L-k, so the average for all 
PNe would be much lower). Considering the relatively 
large uncertainties, the above values show that, for a 
given central star luminosity, the IR-bright selection cri- 
terion would not preferentially select one type of grains 
above others (at least in broad terms). 

At a lower flux limit of 0.5 Jy, and under the as- 
sumption that 20% of the total luminosity is radiated 
in the 12/im IRAS band, a good portion of the galactic 
disk should be sampled: the maximum distance at which 
a PN with warm dust would be detected is 20kpc, for 
L*= 10000 L0. It is thus apparent that the complete- 
ness and homogeneity of the sample discussed here are 
dominated by the selection effects in the PN catalogues 
themselves. As PNe are, for the most part, discovered 
through optical surveys, the distances are unlikely to 
be much in excess of 3 kpc, especially towards the inner 
Galaxy. 



7 THE GALACTIC DISK DISTRIBUTION 
OF PN DUST EMISSION FEATURES 

The spatial distribution of PNe dust types can now be 
constructed using Table ^ and is shown in Figure 6 
in the case of the Zhang (1995) distance scale. Table 6 
lists the properties of the distribution. The total number 
of objects is Ntot=54, consisting of 33 for R < Ro, 
and 21 for R > R^ {Ro = 8.5 kpc, Kerr & Lynden-Bell 
1986). The decrease in the relative proportion of Silicate 
PNe, which was hinted at in Section ^, is confirmed to 
a somewhat higher degree of significance: the fraction 
of Silicate PNe decreases from 0.27±0.08 for R < Ro to 
0.14±0.08 for R > Ro. There is a concentration of UIR 
PNe towards z=0, confirming that they are related to 
higher progenitor masses. On the other hand SiC and 
silicate nebulae have a similar spread in height above 
the galactic plane. However, the vertical distribution of 
PN dust types is rather homogeneous when compared 
to that obtained in the Peimbert types. 

Table n summarises the properties of the PN dis- 
tribution obtained with Djbas distances. It may seem 
surprising that SiC nebulae are at greater distances on 
average, but it is compatible with finding most SiC neb- 
ulae outside the solar circle, where PN catalogues are 
less affected by interstellar extinction. The proportion 
of silicate PNe is confirmed to decrease with i? at a 
higher degree of significance: the fraction of silicate PNe 
decreases from 0.32 ± 0.09 for R < Ro to 0.14±0.06 for 
R > Ro- 

As is apparent from Figure |6|, 'weak' nebulae are 
closer on the Zhang (1995) scale. It was mentioned that 



Table 6: The properties of the distribution of PN dust 
types, based on the distances from Zhang (1995). The 
errors quoted correspond to one standard deviation. 
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PNe with no warm dust are equally distributed among 
Peimbert types, and that their gas phase C/O ratios 
reflect the proportions for the whole sample. They also 
have higher 6 cm fluxes, which together with the selec- 
tion criteria of compact angular size, explains why they 
are on average closer. It is thus very likely that PNe 
with weak continuum correspond to later evolutionary 
stages, and are not a transition stage where C/0~1. 
Although DjuAS distances are not applicable to 'weak' 
PNe, whose optical thickness is uncertain, it is interest- 
ing to note that the average Djuas distance to 'weak' 
PNe is 9.0 kpc, greater than for any other type of PNe. 
This could be interpreted as a lower average luminosity 
(as expected for more evolved PNe on the white dwarf 
cooling track), or that for 'weak' PNe the far-IR flux 
FiRAS is not a good approximation to the bolometric 
flux. 

There is a peculiar asymmetry in the face-on dis- 
tribution of PNe of Figure o. The sector of the galactic 
disk with southern galactic longitudes (the third and 
fourth quadrants) is underpopulated. This is an effect 
due simply to the incompleteness of the catalogues. The 
same asymmetry can be seen in the face-on map of 
Durand et al. (1998), with a larger number of PNe. 
Warm dust PNe with reliable IRAS fluxes all gather 
in very tight IRAS colour-colour boxes (in particular 
log(F(100/im)/F(60Aim))< 0, log(F(25/(m)/F(12/im))> 
0) which allows selecting all warm-dust PN candidates 
from the IRAS PSC. We found 331 IRAS point sources 
with colours of warm dust PNe, whose galactic longi- 
tude/latitude distribution is uniform from northern to 
southern longitudes. Also, the Carina spiral arm be- 
tween the third and fourth quadrant is viewed tangen- 
tially from the sun, thus increasing the interstellar ex- 
tinction and limiting the PN discovery rate. 



8 CONCLUSIONS 

The total number of PNe with 8-13/im spectra has been 
increased to 74 with the inclusion of 24 new objects. 
The sample consists of compact and IR-bright galactic 
disk PNe listed in the Strasbourg-ESO catalogue. 54 
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Table 7: Same as table pi but with Djjias distances. 
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Figure 6: The galactic disk distribution of warm-dust 
PNe on a face-on view of the disk, with distances from 
Zhang (1995). Ticks are every 5 kpc, and the solar circle 
is shown in dashed line. 



PNe have clearly identified warm dust emission features, 
which are placed into three groups (see Table ||): 12 PNe 
show silicate emission, 16 show SiC, and 26 show the 
UIR bands. The remainder have 8-13/im spectra domi- 
nated by emission lines, and correspond to later evolu- 
tionary stages. Thus 22±6% of the PNe with warm dust 
emission have 0-rich grains. We have used this sample 
for an initial study of the PN dust emission features in 
the galactic context. 

A comparison of the PNe dust types with the gas 
phase C/0 ratio shows a good correspondence: Sili- 
cate nebulae have C/0<1, SiC nebulae are found with 
C/0>1, while PNe that show the UIR bands often have 
C/0>>1. We thus confirm that the dust emission fea- 
tures represent an alternative to the plasma diagnostic 
for measuring the C/0 chemical balance in PNe. Neb- 
ulae that show the UIR emission bands also have the 
highest N/0 gas phase ratio. Silicate nebulae are found 
either with high N/0 ratios, or no nitrogen enrichment 
at all. On the other hand SiC nebulae are more uni- 
formly distributed in N/0 ratios. Thus, on a relative 
mass scale, PNe with emission from the UIR bands cor- 
respond to higher progenitor masses, and those with 
SiC to intermediate masses. Silicates are found mainly 
for low mass progenitors, but also for the most massive 
ones. The dust emission features thus provide comple- 
mentary information on the progenitors masses to the 
Peimbert types. 

The adoption of statistical distances showed that 
the sample is large enough to show stratification in 
Peimbert types. We find a link between objects with 
UIR band emission and higher progenitor masses, as 
indicated by their concentration towards the galactic 
plane, obtained from their sky distribution and through 
the use of two independent PN distance scales. There 
is a trend for a decreasing proportion of 0-rich PNe 
with galactocentric radius, confirmed by both distance 
scales, from 30±9% inside the solar circle, to 14±7% 
outside. This trend reflects the variations in the M/C 
star ratio from Thronson et al. (1987) and Jura et al. 
(1989). 

We also showed that the /i?^ 5" fluxes are a good rep- 
resentation of the bolometric flux for PNe with warm- 
dust emission (Section H). The requirement F(12/xm) > 
0.5 Jy should probe a good portion of the galactic disk, 
and the dominant selection effects are rooted in the PNe 
catalogues. 

Although most known IR bright and compact PNe 
were included in this study, further observations are re- 
quired to improve the statistics. Large aperture tele- 
scopes and mid-IR array detectors will be much more 
sensitive for the detection of the warm dust emission 
features in these compact objects, and could allow a 
more accurate analysis. 



12 S. Casassus et al. 



ACKNOWLEDGMENTS 

We are grateful to the referee for interesting comments, 
and to PATT for the time allocation on UKIRT, which 
is operated by the Joint Astronomy Centre on behalf 
of the Particle Physics and Astronomy Research Coun- 
cil, and PATT and ATAC for allocations on the A AT. 
S.C. acknowledges support from Fundacion Andes and 
PPARC through a Gemini studentship. 



REFERENCES 

Acker A., Ochsenbein F., Stcnholm B., Tylenda R., Marcout J., 
Schohn C, 1992, 'The Strasbourg-ESO Catalogue of Galac- 
tic Planetary Nebulae', published by the European Southern 
Observatory 
Aitkcn D.K., Barlow M.J., Roche P.F., Spenser P.M., 1980, MN- 

RAS, 192, 679 
Aitken D.K., Roche P.F., 1982, MNRAS, 200, 217 
Aitken D.K., Roche P.F., Spenser P.M., Jones B., 1979, ApJ, 233, 

925 
Aller L.H., Czyzak S.J., 1983, ApJSS, 51, 211 
Aller L.H., Hyung S., 1995, MNRAS, 276, 1101 
Aller L.H, Hyung S., Feibelman W.A., 1996, PASP, 108, 488 
Aller L.H., Keyes CD., 1987, ApJSS, 65, 405 
Aller L.H., Ross J.E., O'Mara B.J., Keyes CD., 1981, MNRAS, 

197,95 
Andersen A.C, Jager C, Mutschke H., Braatz A., Clement D., 
Henning Th., j0rgensen U.G., Ott U., 1999, A&A, 343, 933 
Barlow M.J., 1983, lAU symp 103 "Planetary Nebulae" ed D.R. 

Flower pl05 
Blocker T., 1995, A&A, 299, 755 

Cahn J.H., Kaler J.B., StangheUini L., 1992, A&AS, 94, 399 
Casassus S., Roche P.F., 2000, submitted (paper II) 
Cohen M., AUamandola L., Tielens A.G.G.M., Bregman J., Simp- 
son J.P., Witteborn F.C, Wooden D., Rank D.M., 1986, ApJ, 
302, 737 
Cohen M., Jones B.F., 1987, ApJ 321, L151 

Cohen M., Tielens A.G.G.M., Bregman J., Witteborn F.C, Rank 
D.M., AUamandola L.J., Wooden D.H., dc Muizon M., 1989, 
ApJ, 341, 246 
Costa R.D.D., Chiapini C, Maciel W.J., de Freitas Pacheco J. A., 

1996, AASS, 116, 249 
Crowther P.A., De Marco O., Barlow M.J., 1998, MNRAS 296, 

367 
Cuisinier F., Acker A., Koppen J., 1996, A&A 307, 215 
Day K.L., 1979, ApJ, 234, 158 
Day K.L., 1981, ApJ, 246, 110 
Day K.L., Donn B., 1978, ApJ, 222, L45 

De Marco O., Barlow M.J., Storey P.J., 1997, MNRAS, 292, 86 
Duley W.W, Williams D.A., 1981, MNRAS 196, 269 
Durand S., Acker A., Zijlstra A., 1998, AAS, 132, 13 
Faundez-Abans M., Maciel W.J., 1987, A&A, 183, 324 
Feibelman W.A., Aller L.H, Keyes CD., Czyzak S.J., 1985, Proc. 

Nat. Acad. Sci., 82, 2202 
Feibelman W.A., Hyung S., Aller L.H., 1994, ApJ, 426, 653 
Forrest W.J., Gillett F.C, Stein W.A., 1975, ApJ, 195, 423 
Gillett F.C, Forrest W.J., Merrill K.M., 1973, ApJ, 183, 87 
Gillett E.G., Jones T.W., Merrill K.M., Stein A.W., 1975, ApJ, 

45, 77 
Gilman R.C, 1969, MNRAS, 302, 305 

Gomez Y., Rodriguez L.F., Moran J.M., 1993, ApJ, 416, 620 
Grevesse N., Anders E., 1989, 'Cosmic Abundances of Matter', 
AIP Conference Proceedings, 183, 1, Ed. C.J. Waddington, 
(New York: AIP) 
Groenewegen M.A.T, dejong T., 1993, A&A, 267, 410 (GJ93) 



Hajian, A.R., Terzian Y., Bignell C, 1995, AJ, 109, 2600 

Hajian, A.R., Terzian Y., 1996, PASP, 108, 258 

Hyung S., 1994, ApJS 90, 119 

Hyung S., Aller L.H., 1995, NNRAS, 273, 257 

Hyung S., Aller L.H., 1996, MNRAS, 278, 551 

Hyung S., Aller L.H., 1997, MNRAS, 292, 71 

Hyung S., Aller L.H., Feibelman W.A., 1994a, MNRAS, 269, 975 

Hyung S., Aller L.H., Feibelman W.A., 1994b, ApJSS, 94, 465 

Hyung S., Aller L.H., Feibelman W.A., 1994c, PASP, 106, 745 

Hyung S., Aller L.H., Feibelman W.A., 1997, ApJSS, 108, 503 

Hyung S., Keyes CD., Aller L.H., 1995, MNRAS, 272, 49 

Jura M., Joyce R.R., Kleinmann S.C, 1989, ApJ, 336, 924 

Kerr F.J., Lynden-Bell D., 1986, MNRAS, 221, 1023 

Keyes CD., Aller L.H., Feibelman W.A., 1990, PASP, 102, 59 

Kingsburgh R.L., Barlow M.J., 1994, MNRAS, 271, 257 

Koppen J., Acker A., Stcnholm B., 1991, A&A, 248, 197 

Kwok S., 1994, PASP, 106, 344 

Kwok S., Hrivnak B.J., Milone E.F., 1986, ApJ, 303, 451 

Leger A., Puget J.L., 1984, A&A, 137, 5 

Leisy P., Dennefeld M., 1996, A&ASS, 116, 95 

Likkcl L., Morris M., Kastner J.H., Forveille T., 1994, A&A, 282, 
190 

Lutz J.H., 1984, ApJ 279, 714 

Maciel W.J., Dutra CM., 1992, A&A 262, 271 

Mendez R.H., Kudritzki R.P., Herrero A., 1992, A&A, 260, 329 

Mishurov, Y.N., Zenina I.A., 1999, A&A 341, 81 

Paczynski B., 1970, Acta Astron., 20, 47 

Peimbert M., 1978, in Terzian Y., ed., Proc. lAU Symp 76, 'Plan- 
etary Nebulae'. Reidel, Dordrecht, p. 233 

Peimbert M., 1992, in 'Observational Astrophysics', Ed. R.E. 
White, lOP Publishing Ltd. 

Perinotto M., 1991, ApJSS 76, 687 

Roche P.F., 1989, lAU symp 131 "Planetary Nebulae" cd S.Torres 
Peimbert, pll7 

Roche P.F., Allen D.A., Aitken D.K., 1983a, MNRAS, 204, 1009 

Roche P.F., Aitken D.K., 1983, MNRAS, 203, 9 

Roche P.F., Aitken D.K., 1986, MNRAS, 221, 63 

Roche P.F., Aitken D.K., Smith C.H., 1991, MNRAS, 252, 282 

Roche P.F., Aitken D.K., Whitmore B., 1983b, MNRAS, 204, 
1017 

Roche P.F., Lucas P.W., Hoare M.G., Aitken D.K., Smith C.H., 
1996, MNRAS 280, 924 

Russell R.R., Soifer B.T., Willncr S.P., 1978, ApJ, 191, 141 

Salpeter E.E., 1955, ApJ 121, 161 

Smith, C.H., Aitken, D.K., Roche P.F., Wright, CM., 1995, MN- 
RAS, 277, 259 

Sylvester R.J., Skinner C.J., Barlow M.J., 1998, MNRAS, 301, 
1083 

Thronson H.A., Latter W.B, Black J.H., Bally J., Hacking P., 
1987, ApJ, 322, 770 

Van de Steene G.C, Zijlstra A.A., 1994, AASS, 108, 485 

Wagenhuber J., Groenewegen M.A.T. , 1999, A&A, 340, 183 

Willner S.P., Jones B., Puetter R.C, Russell R.W., Soifer B.T, 
1979, ApJ, 234, 496 

Zhang C.Y., 1993, ApJ, 410, 239 

Zhang C.Y., 1995, ApJSS, 98, 659 

Zhang C.Y., Kwok S., 1993, ApJSS, 88, 137 

Zuckerman B., Aller L.H., 1986, ApJ, 301, 772 

This paper has been produced using the Royal Astro- 
nomical Society/Blackwell Science MJrjX style file. 



